Orbital fluctuations in the RVO3 perovskites 

Andrzej M. Oles and Peter Horsch 



Abstract The properties of Mott insulators with orbital degrees of freedom are de- 
scribed by spin-orbital superexchange models, which provide a theoretical frame- 
work for understanding their magnetic and optical properties. We introduce such a 
model derived for (xy) 1 (yz/zx) 1 configuration of V 3+ ions in the perovskites, 
/?=Lu,Yb,- • - ,La, and demonstrate that {yz,zx} orbital fluctuations along the c axis 
are responsible for the huge magnetic and optical anisotropics observed in the al- 
most perfectly cubic compound LaV03. We argue that the GdFe03 distortion and 
the large difference in entropy of C-AF and G-AF phases is responsible for the sec- 
ond magnetic transition observed at T^2 in YVO3. Next we address the variation of 
orbital and magnetic transition temperature, 7bo and Tjyi , in the R VO3 perovskites, 
after extending the spin-orbital model by the crystal-field and the orbital interactions 
which arise from the GdFeC>3 and Jahn-Teller distortions of the VO6 octahedra. We 
further find that the orthorhombic distortion which increases from LaVC>3 to LUVO3 
plays a crucial role by controlling the orbital fluctuations, and via the modified or- 
bital correlations influences the onset of both magnetic and orbital order. 



1 Orbital degrees of freedom in strongly correlated systems 

Orbital degrees of freedom play a key role for many intriguing phenomena in 
strongly correlated transition metal oxides, such as the colossal magnetoresistance 
in the manganites or the effective reduction of dimensionality in KCUF3 QJ. Before 
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addressing complex phenomena in doped Mott insulators, it is necessary to describe 
first the undoped materials, such as LaMnC>3 or LaVC>3. These two systems are 
canonical examples of correlated insulators with coexisting magnetic and orbital or- 
der JT|. In both cases large local Coulomb interaction U suppresses charge fluctua- 
tions, leading to low-energy effective Hamiltonians with superexchange interactions 
which stabilize antiferromagnetic (AF) spin order at low temperature J2][3|. How- 
ever, the AF order is different in both cases: ferromagnetic (FM) planes are coupled 
by AF interactions in the A -type AF phase of LaMn03, while FM chains along the c 
cubic axis are coupled by AF interactions in the ab planes in the C-type AF (C-AF) 
phase of LaV03. The superexchange Hamiltonians which describe both systems are 
just examples for the spin-orbital physics J4), where orbital (pseudospin) operators 
contribute explicitly to the structure of the superexchange interactions — their actual 
form depends on the number of 3d electrons (holes) at transition metal ions which 
determines the value of spin S, and on the type of active orbital degrees of freedom, 
e g or tig. In simple terms, the magnetic structure is determined by the pattern of 
occupied and empty orbitals, and the associated rules are known as Goodenough- 
Kanamori rules (GKR). The central focus of this overview are f? g orbital degenerate 
systems, where quantum fluctuations of orbitals play a central role for the electronic 
properties |3]|5]|6) and modify the predictions of the GKR. 

In the last two decades several new concepts were developed in the field of or- 
bital physics |fl~). The best known spin-orbital superexchange Hamiltonian is the 
Kugel-Khomskii model [7|, which describes the e g orbital {x 2 — y 2 ,3z 2 — r 2 } de- 
grees of freedom coupled to S = 1/2 spins at Cu 2+ (d 9 ) ions in KCUF3. The spins 
interact by either FM and AF exchange interactions, depending on the type of occu- 
pied and empty orbitals on two neighboring ions. It has been found that enhanced 
quantum fluctuations due to orbital degrees of freedom, which contribute to joint 
spin-orbital dynamics, may destabilize long-range magnetic order near the quantum 
critical point of the Kugel-Khomskii model |8). The orbital part of the superex- 
change is thereby intrinsically frustrated even on geometrically non-frustrated lat- 
tices, as in the perovskite lattice H|9), which is a second important concept in the 
field of orbital physics. Finally, although spin and orbital operators commute, there 
are situations where joint spin-orbital dynamics plays a crucial role, and spin and 
orbital operators cannot be separated from each other. This situation is called spin- 
orbital entanglement ifTol . and its best example are the entangled SU(4) singlets in 
the one-dimensional (ID) SU(4) model [JJJ. There is no doubt that these recent de- 
velopments in the orbital physics provide many challenges both for the experimental 
studies and for the theoretical understanding of the experimental consequences of 
the spin-orbital superexchange. 

Let us consider first the orbital part of the superexchange. Its intrinsic frustra- 
tion results from the directional nature of orbital pseudospin interactions flUED — 
they imply that the pair of orbitals which would minimize the energy depends on 
the direction of a bond (ij) in a cubic (peovskite) lattice. In case of e g orbitals 
the superexchange interactions are Ising-like as only one orbital flavor allows for 
electron hopping t and the electron exchange process does not occur. They favor a 
pair of orthogonal orbitals on both sites of the considered bond (fl~2), for instance 
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\z) ~ (3z 2 — r 2 )/ Vo and |x) ~ x 2 — y 2 orbital for a bond along the c axis. When the 
two above orbital states are represented as components of x = 1/2 pseudospin, this 
configuration gives the energy of — where J is the superexchange constant. Un- 
like in the ID model fPHl . such an optimal orbital configuration cannot be realized 
simultaneously on all the bonds in a two-dimensional (2D) or three-dimensional 
(3D) system. Thus, in contrast to spin systems, the tendency towards orbital disor- 
dered state (orbital liquid) is enhanced with increasing system dimension II14II15I . 

The essence of orbital frustration is captured by the 2D compass model, origi- 
nally developed as a model for Mott insulators Q. Intersite interactions in the com- 
pass model are descibed by products t"t™ of single pseudospin components, 

*f = 5°?, ^ = 1(7?, *1 = \°!- (1) 

for a bond (ij) | 7, where a — x,y,z, rather than by a pseudospin scalar product 
Tj • Xj. For instance, in the 2D case of a single ab plane, the compass model Ifl6l , 

H 2D =J X £ T?X*+J Z £ TfTj, (2) 

<<7)ll« <y>H6 

describes the competition between xfXj and xfxj interactions for the bonds along a 
and b axis, respectively. This competition of pseudospin interactions along differ- 
ent directions results in intersite correlations similar to those in the anisotropic XY 
model, and generates a quantum critical point at J x = with high degeneracy of the 
ground state ifTTl . So, despite certain similarities of the compass model to ordinary 
models used in quantum magnetism, an ordered phase with finite magnetization is 
absent. It is interesting to note that a similar quantum phase transition exists also in 
the ID chain compass model IT8l (N' = N/2 is the number of unit cells): 

N' 

Hw = £ {/tifi-i^/ + 4:44+1 } • (3) 
1=1 

Recently this ID compass model was solved exactly in the whole range of {J X ,J Z } 
parameters [ 18] by mapping to the exactly solvable (quantum) Ising model in trans- 
verse field. It provides a beautiful example of a first order quantum phase transition 
between two phases with large (4-14) or (44+i) correlations, and discontinu- 
ous changes of intersite correlation functions. 

In realistic spin-orbital superexchange models transitions between different or- 
dered or disordered orbital states are accompanied by magnetic transitions. This 
field is very rich, and several problems remain unsolved as simple mean-field (MF) 
approaches do not suffice in general, even for the systems with perovskite lat- 
tices |@). In this chapter we shall address the physical properties of the RVO3 per- 
ovskites (7?=Lu,Yb,- • - ,La), where not only the above intrinsic frustration of the or- 
bital superexchange, but also the structure of the spin-orbital superexchange arising 
from multiplet splittings due to Hund's exchange plays a role and determines the ob- 
served physical properties at finite temperature. Moreover, we shall see that the cou- 
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Fig. 1 Magnon dispersion 
relation obtained by neutron 
scattering for the C-AF phase 
of YVO3 at T = 85 K. The 
lines are interpolation be- 
tween the experimental points 
(squares with error bars) along 
two high symmetry directions 
in the Brillouin zone. Image 
courtesy of Clemens Ulrich. 
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pling of the orbitals to the lattice, i.e., via Jahn-Teller (JT) coupling, GdFeOg-like 
and orthorhombic distortion, are important control parameters. First we analyze the 
structure of the spin-orbital superexchange in section 2 and show its consequences 
for the magnetic and optical properties of strongly correlated transition metal com- 
pounds. Here we also address the entanglement of spin and orbital variables which 
is ignored in the MF decoupling, and we point out that it fails in certain situations. 

The coupling between the orbital and spin variables is capable of generating qual- 
itatively new phenomena which do not occur in the absence of orbital interactions, 
such as anisotropic magnetic interactions, and novel quantum phenomena at finite 
temperature, discusseed on the example of LaVC>3 in section 3. One of such novel 
and puzzling phenomena is the magnetic phase transition between two different 
types of AF order observed in YVO3 — this compound has G-type AF order (stag- 
gered in all three directions, called below G-AF phase) at low temperature T < Tmi, 
while the magnetic order changes in the first order magnetic transition at 7)v2 = 77 
K to C-AF phase which remains stable up to Ty\ ~ 1 16 K. The latter C-AF phase 
has rather exotic magnetic properties, and the magnon spectra show dimerization 
of the FM interactions along the c axis 1191 . see figure Q] In fact, the G-AF phase 
occurs in systems with large GdFe03-like distortion [20|. In Ref. (3) an orbital in- 
teraction favoring C-type alternating orbital (C-AO) order was invoked to explain 
the G-AF phase. We also address this problem in section 3 and present arguments 
that at higher T > C-AF phase reappears is due to its higher entropy J3|. 

In section 4 we address the experimental phase diagram of the perovskites. 
It is quite different from the (also puzzling) phase diagram of the /?Mn03 per- 
ovskites, where the orbital order (OO) appears first at Too upon lowering the temper- 
ature, and spin order follows 12"T1 at the Neel temperature, T^i -C 7bo- In contrast, in 
the RVO3 perovskites the two transitions appear at similar temperature 11221 . For in- 
stance, in LaV03 they occur even almost simultaneously i.e., T^i ~ Too- However, 
they become separated from each other in the R VO3 systems with smaller ionic radii 
of R ions — whereas 7]yi gets reduced for decreasing r^, Too exhibits a nonmono- 
tonic dependence on r# ll20l . A short summary is presented in section 5. We also 
point out a few unsolved problems of currect interest in the field of orbital physics. 
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2 Spin-orbital superexchange and entanglement 

Spin-orbital models derived for real systems are rather complex and the orbital part 
of the superexchange is typically much more complicated than the compass model 
(f2]i of section 1 10] . The main feature of these superexchange models is that not 
only the orbital interactions are directional and frustrated, but also spin correlations 
may influence orbital interactions and vice versa. This is best seen in the Kugel- 
Khomskii (d 9 ) model |7j, where the G-AF and A-AF order compete with each other, 
and the long-range order is destabilized by quantum fluctuations in the vicinity of 
the quantum critical point (Jh,E z ) = (0, 0) |23l . Here Jh is the local exchange (see 
below), and E z is the splitting of two e g orbitals. Although this model is a possible 
realization of disordered spin-orbital liquid, its phase diagram remains unexplored 
beyond the MF approach and simple valence-bond wave functions of Ref. 1 8 1 — it 
remains one of the challenges in this field. 

In this chapter we consider the superexchange derived for an (idealized) per- 
ovskite structure of RVOj,, with V 3+ ions occupying the cubic lattice. The kinetic 
energy is given by: 

ff/ = -'E L ( d lac d jaa+d) aa d iac ), (4) 

<<7}||y a(r)>°- 

where dj aa is electron creation operator for an electron with spin a =t, 4- m orbital 
a at site ;. The summation runs over the bonds (ij) \\y along three cubic axes, 7 = 
a,b,c, with the hopping elements t between active orbitals. They originate from 
two subsequent hopping processes via the intermediate 2p K oxygen orbital along 
each V-O-V bond. Its value can in principle be derived from the charge-transfer 
model ll24ll . and one expects t = t 2 d /A ~ 0.2 eV J3). Only two out of three ti g 
orbitals, labelled by oc(y), are active along each bond (ij) and contribute to the 
kinetic energy (O, while the third orbital lies in the plane perpendicular to the 7 axis 
and the hopping via the intermediate oxygen 2p K oxygen is forbidden by symmetry. 
This motivates a convenient notation used below, 

|a) = |yz), \b) = \xz), \c) = \xy), (5) 

where the orbital inactive along a cubic direction 7 is labelled by its index as \y). 

The superexchange model for the RWO3 perovskites arises from virtual charge 
excitations between V 3+ ions in the high-spin S=l state. The number of d electrons 
is 2 at each V 3+ ion (d 2 configuration), and the superexchange is derived from all 
possible virtual dfd 2 f± dfdj excitation processes (for more details see Ref. l25l ). 
It is parametrized by the superexchange constant J and Hund's parameter 77, 

j 4f2 J " 

J = — , n = — , (6) 

U U 

where U is the intraorbital Coulomb interaction and Jh is Hund's exchange between 
t% R electrons. Here we use the usual convention and write the local Coulomb inter- 
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actions between 3d electrons at V 3+ ions by limiting ourselves to intraorbital and 
two-orbital interaction elements J4): 

#int = ^X«iat%4+ (u- ~ J h) Y, n ian^-2J H Y Sia-Sifi 
ia i,a<p i.a<p 

+ j h L ( d U d U d iPi d ip^+ d m d lpi d i^ iat ). (7) 

i,a</3 

When only one type of orbitals is party occupied (as in the present case of the RVO3 
perovskites or in KCUF3), the two parameters {U,Jh} are sufficient to describe these 
interactions in Eq. (0: (;) the intraorbital Coulomb element U and (;;) the interor- 
bital (Hund's) exchange element Jh, where {A,B,C} are the Racah parameters l26l . 
In such cases the above expression is exact; in other cases when both e g and t% g elec- 
trons contribute to charge excitations (as for instance in the /?Mn03 perovskites), 
Eq. (0 is only an approximation — the anisotropy on the interorbital interaction 
elements has to be then included to reproduce accurately the multiplet spectra of the 
transition metal ions |26 |. The intraorbital interaction is C/ = A +4-B + 3C, while Jh 
depends on orbital type — for ti„ electrons one finds [4 26 1 Jh = 3B + C. 

The perturbative treatment of intersite charge excitations dfd 2 =± dfdj in the 
regime of t <C U leads for the RVO3 perovskites (and in each similar case 0]) to the 
spin-orbital superexchange model: 

m= E H^(ij)=J £ {{SrS J+ S 2 )j)p+K^}. (8) 

<y)llr (y)llr 

The spin interactions « S, • S ; - obey the SU(2) symmetry. In contrast, the orbital 

interaction operators jjj and iDv involve directional (here t2 g ) orbitals on each 
individual bond (ij) \\ 7, so they have a lower (cubic) symmetry. The above form of 
the spin-orbital interactions is general and the spin value S depend on the electronic 
configuration d" of the involved transition metal ions (here n = 2 and 5=1). For 

convenience, we introduced also a constant S 2 in the spin part, so for the classical 

Hy) 

Neel order the first term <x J>!' vanishes. 

In the R VO3 perovskites one finds the orbital operators l25l : 

Jl? = \{(l+2nn)(r i -r j + \n i n J j 

- 77 r 3 f T; x Xj + \n in ^\ -\nn («i + nj) J , 



(9) 



►(7) 



1 



K>!> =\X]n[ Tj • Tj + -n t n } ) + T]r 3 [ T; x x 



(7) 



-{l + nn ){ ni + nj)\ , (10) 
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where the scalar product (t, • Tjp 7 ' and the cross-product, 

(x l XX j ^ = ^{x+X+ + X7Tj)+XtT j , (11) 

involve orbital (pseudospin) operators corresponding to two active t2 g orbitals along 
the 7 axis, with T; = {T+, t^t^}, and 

if = ±(n Uyz -n iiZX ). (12) 

They follow from the structure of local Coulomb interaction ©. The latter term (fTTT i 
leads to the nonconservation of total pseudospin quantum number. Density operators 
nf" 1 in Eqs. © and ( fTOb stand for the number of d electrons in active orbitals for the 
considered bond (ij), e.g. nf' = n\ a +«,;,. The coefficients, 

ri = i^W r3 = r+W (13) 

follow from the energies of dfd 2 ^ dfdj excitations in the units of U : (;) n repre- 
sents the high-spin 4 A2 excitation of energy (U — 3Jh), while the low-spin excita- 
tions are given by (z'z) r^ = \ originating from the low-spin 2 T\ and 2 E excitations 
of energy U, and {Hi) rj, represents the low-spin 2 72 states of energy (U + 27//). 

Magnetic order observed in Mott insulators is usually understood in terms of 
the GKR which are based on the MF picture and ignore entangled quantum states. 
These rules state that the pattern of occupied orbitals determines the spin structure. 
For example, for 180° bonds (e.g. Mn-O-Mn bonds in LaMnOs) there are two 
key rules: (/) if two partially occupied 3d orbitals point towards each other, the 
interaction is AF, however, (//) if an occupied orbital on one site has a large overlap 
with an empty orbital on the other site of a bond (ij), the interaction is weak and 
FM due to finite Hund's exchange. This means that spin order and orbital order are 
complementary — ferro-like (uniform) orbital (FO) order supports AF spin order, 
while AO order supports FM spin order. Indeed, these celebrated rules are well 
followed in LaMnC<3 11271 and in KCUF3 11281 . where strong JT effect stabilizes the 
orbital order and suppresses the orbital fluctuations. The AO order is here robust in 
the FM ab planes, while the orbitals obey the FO order along the c axis, supporting 
the AF coupling and leading to the A-AF phase for both systems. In such cases the 
GKR directly apply. Therefore, one may disentangle the spin and orbital operators, 
and it has been shown that this procedure is sufficient to explain both the magnetic 
(2) and optical J4] properties of LaMn03. 

As another prominent example of the Goodenough-Kanamori complementarity 
we would like to mention the AF phases realized in YVO3 lfT9l . which are the sub- 
ject of intense research in recent years. A priori, the orbital interactions between 
V 3+ ions in d 2 configuration obey the cubic symmetry, if the t2 g orbitals are ran- 
domly occupied. However, the symmetry breaking at the structural transition where 
the symmetry is reduced from cubic to orthorhombic, which persists in the magnetic 
phases, suggests that the electronic configuration is different. Indeed, the GdFe03 
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distortions in the RMOj, structure break the symmetry in the orbital space, and both 
the electronic structure calculations ||29l and the analysis using the point charge 
model [30 1 indicate that the electronic configuration (xy) 1 (yz,zx) 1 is induced at ev- 
ery site, i.e., 

«, c = l, n ia +n ib = \. (14) 

The partly filled {a,b} orbitals are both active along the c axis, and may lead either 
to FO or to AO order. Indeed, depending on this orbital pattern, the magnetic cor- 
relations are there either AF or FM, explaining the origin of the two observed types 
of AF order: (;) the C-AF phase, and (ii) the G-AF phase. However, the situation is 
more subtle as both orbitals in the orbital doublet { \yz) , \xz)} = \b)} at each site 
i are active on the bonds along the c axis. This demonstrates an important difference 
between the e g (with one electron or one hole in active e ? orbital at each site |fl3l ) 
and a t2 g system, such as RVOt, perovskite vanadates, where electrons occupying 
two active t2 g orbitals may fluctuate and form an orbital singlet (3). The cubic sym- 
metry is thus broken as both orbital flavors are active only along the c axis, and 
the bonds in the ab planes and along the c axis are nonequivalent. Consequently, 
superexchange orbital operators © and dTob take different forms along these two 
distinct directions, 

j\f = \ {(1 +2ijn) U ■ xj + T?r 3 U x tj + - r]n\ , (15) 

l if = j 7 ?'"! ( T <-Ty + ^+77r 3 (xiXXj + ^j -1(1 + 77^1 , (16) 

4? = \{{l-n^){\+n lb n ]b )- n {n lb -n jb ) 2 } , (17) 

4 0) = \n{nn+ ri ){\+n lb n ]b ). (18) 

The general form of spin-orbital superexchange model (HJ suggests that the above 
symmetry breaking leads indeed to an effective spin model with broken symmetry 
between magnetic interactions along different cubic axes. By averaging over the or- 
bital operators one finds indeed different effective magnetic exchange interactions, 
J c along the c axis and J ab within the ab planes: 

Jc = {S\f), J ab = (fi?)- d9) 

The interactions in the ab planes could in principle still take two different values in 
case of finite lattice strain discussed below, making both {a,b} axes inequivalent, 
but here we want just to point out the symmetry breaking between the c axis and 
the ab planes, which follows from the density distribution ( [Pfl ) and explains the 
nonequivalence of spin interactions in the C-AF phase of the RVO3 perovskites Q. 

Apart from the superexchange there are in general also interactions due to the 
couplings to the lattice that control the orbitals. In the cubic vanadates these inter- 
actions are expected to be weak, but nevertheless they influence significantly the 
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spin-orbital fluctuations and decide about the observed properties in the RVO3 fam- 
ily. We write the orbital interactions, « tfXj, induced by the GdFe03 distortions and 
by the JT distortions of the lattice using two parameters, V a b and V c , 

^v=V ab £ TfTj-V c £ Tfrj. (20) 

(>j)\\c (ij)\\c 

The orbital interaction along the c axis V c plays here a crucial role and allows one 
to switch between the two types of magnetic order, C-AF and G-AF phase OTIl . 
stabilizing simultaneously either G-AO or C-AO order. 

However, the description in terms of the GKR does not suffice and the ground 
state of spin-orbital model for the RVO3 perovskites, which consists of the superex- 
change and the effective orbital interactions, 

■ffl$ z = Mj + My . (21) 

may also be entangled due to the quantum coupling between spin 5=1 and orbital 
1 = 1/2 operators along the c axis, see Eq. (TT~5T >. In constrast, the orbital fluctua- 
tions in the ab planes are quenched due to the occupied c orbitals at each site ( TBI ), 
so spins and orbitals disentangle. Possible entanglement between spin (S; • S ; ) and 
orbital (tj • xf) operators along the bonds (if) \\ c in the RVO^ perovskites, and the 
applicability of the GKR to these systems, may be investigated by evaluating in- 
tersite spin and orbital correlations (to make these two functions comparable, we 
renormalized the spin correlations by the factor 3), 

Sy = i<S,.Sy), (22) 
Tij = (t, • T j) , (23) 

and comparing them with each other. A key quantity that measures spin-orbital en- 
tanglement is the composite correlation function ifTUl . 

Q/=i{((SrSy)(T i .Ty)>-(S,.S / )<T i .Ty>} . (24) 

When Cjj = 0, the spin and orbital operators are disentangled and their MF decou- 
pling is exact, while if Cy < — spin and orbital operators are entangled and the 
MF decoupling not justified. 

The numerical results for a ID chain along the c axis described by vanadate spin- 
orbital model (fJTJ are shown in Fig. [2] One finds entangled spin-orbital states with 
all three Sy, Tji and Cjj correlations being negative in the spin-singlet (S = 0) regime 
of fluctuating yz and zx orbitals, obtained for rj < 0.07 [Fig.|2a)]. Therefore, the 
complementary behavior of spin d22l) and orbital d23l correlations is absent in this 
regime of parameters and the GKR are violated. In addition, composite spin-orbital 
correlations d24b are here finite (Cy < 0), so spin and orbital variables are entangled, 
and the MF factorization of spin-orbital operators fails. In a similar d l model for the 
perovskite titanates (with S = 1/2) one finds even somewhat stronger spin-orbital 
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Fig. 2 Evolution of intesite correlations and exchange constants along the c axis obtained by exact 
diagonalizaton of spin-orbital model on a chain of N = 4 sites with periodic boundary conditions, 
with ff) and k\f given by Eqs. {B} and {16), for increasing Hund's exchange 77: (a),(b) intersite 
spin Sjj j22t (filled circles), orbital )23b (empty circles), and spin-orbital C, ; j24t (x) correlations; 
(c),(d) the corresponding spin exchange constants 7 I; (19) . In the shaded areas of (c) and (d) the 
spin correlations 5; ; < do not follow the sign of the exchange constant 7, ; < 0, and the classical 
GKR are violated. Parameters: (a),(c) V c = 0, and (b),(d) V c = J. 



entanglement and the regime of 77 with Cy < is broader (i.e., 77 < 0.21) iflOl . At 
the point 77 = one recovers then the SU(4) model with Sy = Tij = Cij = —0.25, 
and the ground state is an entangled SU(4) singlet, involving a linear combination 
of (spin singlet/orbital triplet) and (spin triplet/orbital singlet) states. 

To provide further evidence that the GKR do not apply to spin-orbital model 
d2TT i in the regime of small 77, we compare spin exchange constants Jn ( fT~9b shown 
in Fig. [2c) with spin correlations Sy ( |22] |. see Fig. |2ja). One finds that exchange 
interaction is formally FM (/y < 0) in the orbital-disordered phase in the regime 
of 77 < 0.07, but it is accompanied by AF spin correlations (5y < 0). Therefore 
JijSij > and the ground state energy would be enhanced in an ordered state, 
when calculated in the MF decoupling of spin-orbital operators iflOll . This at first 
instance somewhat surprising result is a consequence of 'dynamical ' nature of ex- 
change constants T>.' which exhibit large fluctuations [10|, measured by the sec- 
ond moment, 8J = {((jjf) 2 ) ~^y} 1//2 - F° r instance, in d 2 model (|2TT > the or- 
bital bond correlations change dynamically from singlet to triplet, resulting in large 
8J= \{l - {ITij+lf} 1 ' 2 ~ 0.247, i.e., 8J > |/y|. 
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Remarkably, finite spin-orbital correlations Q; < and similar violation of the 
GKR are found also at finite orbital interaction d20l i induced by the lattice, V c > 0. 
Representative results obtained for V c = J are shown in Figs.|2jb) and|2jd). At small 
77 FO order is induced, and in this regime the GKR are followed by the AF/FO phase 
(similar to the FM/AO phase at large 77 which also follows the GKR). However, for 
intermediate Hund's exchange rj ~ 0.07 FO order is destabilized and the entangled 
AF/AO phase appears, with similar spin, orbital and composite spin-orbital coor- 
relations as found before at V c = and tj = [Figs. 12a)]. Also in this case FM 
exchange (Jij < 0) coexists with AF spin correlations < 0). Thus we conclude 
that orbital interactions induced by the lattice modify the regime of entangled spin- 
orbital states in the intermediate AF/AO phase which may be moved to more realis- 
tic values of 77 ~ 0.1, and cannot eliminate it completely. In addition, the transition 
between the FO/AF and AO/AF phase is continuous iflOl due to the structure of or- 
bital superexchange which contains terms (fTTT i responsible for non-conservation of 
orbital quantum numbers. 



3 Experimental evidence of orbital fluctuations in LaV03/YVC>3 

Before discussing the exotic magnetic properties and the phase diagram of the R VO3 
perovskites, we will consider the influence of magnetism on the optical spectra of 
LaVOs, starting with a general formulation of the theory. While exchange constants 
may be extracted from the spin-orbital superexchange model dT9l l, it is frequently 
not realized that virtual charge excitations that contribute to the superexchange are 
responsible as well for the optical absorption, thus the superexchange and the optical 
absorption are intimately related to each other via the optical sum rule l33l . This is 
not so surprising as when electrons are almost localized in a Mott insulator, the 
only kinetic energy which is left and decides about the optical spectral weight is 
associated with virtual excitations contributing to superexchange. Therefore, in Mott 
insulators the thermal evolution of optical spectral weight can be deduced from the 

superexchange 041 . In a system with orbital degeneracy the optical spectra consist 

(y) 

of several multiplet transitions, and the kinetic energy K„ (due tod — d excitations) 
associated with each of them can be determined from the superexchange (© using 
the Hellman-Feynman theorem ll32l . 



Note that K„ ' is negative and corresponds to the n'th multiplet state of the transition 
metal ion, created by a charge excitation along a bond (ij) \\ 7. It is obvious that the 
thermal excitation values (■■■) depend sensitively on the magnetic structure, i.e., 
whether spin correlations on a bond (ij) are FM or AF. 

Thus it is natural to decompose the optical sum rule which is usually formulated 
in terms of the total kinetic energy for polarization 7, 




(25) 
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KM=2jZ(Hp\ij)), (26) 

n 

into partial optical sum rules for individual Hubbard subbands 11321 . 

£a^(co)da> = -f & = -x(l#>M) , (27) 

where o„ (©) is the contribution of band n to the optical conductivity for polar- 
ization along the y axis, oq is the distance between transition metal ions, and the 
tight-binding model with nearest neighbor hopping is assumed. Equation ( f27T > pro- 
vides a practical way of calculating the optical spectral weights from spin-orbital 
superexchange models, such as the one derived for the R~VOj, perovskites (|211) , Note 
that the total optical intensity d26l i is of less interest here as it has a much weaker 
temperature dependence and does not allow one for a direct insight into the nature 
of the electronic structure. In addition, it might be also more difficult to resolve from 
experiment. 

In order to apply the above theory to the RVO3 perovskites, we write the su- 
perexchange operator H^\ij) for a bond (ij) | 7, contributing to operator ([8]), 
as a superposition of dfdj ^ dfdj charge excitations to different upper Hubbard 
subbands labelled by n 11321 . 

HM(ij)=£Hp\ij). (28) 

n 

One finds the superexchange terms Hn {ij) for a bond (ij) along the c axis 
Hl c \ij) = ~Jr 1 (2 + S i -S j )(j l 
Hi C \ij) = ~J(l-SvS j )(J i - 
H 3 (c) (y) = -^ 3 (i-S i -S y )(i 

and Hj, ab \ij) for a bond in the ab plane, 

H[ ah) (U) =-i/r 1 (2 + S,-S 7 )Q-TfT?) , (32) 
< ) (';) = -^(l-SrS 7 )(^T^T^-I^), (33) 



- T <- • tj) » 




(29) 






(30) 




• 


(31) 



H <<*)(y) = -- Jn (1 -SrSj) ^ T -,f T -tj + T?t?j . (34) 

These expressions show that the spin correlations along the c axis and within the ab 
planes, 
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Fig. 3 Exchange constants 
J„(, and — J c $19) calculated 
from Eqs. {T7J and H51 in 
the C-AF phase of LaV0 3 
for increasing f] (solid lines). 
Dashed line shows the value 
of — J c obtained for classical 
orbital order <36b according 
to GKR, (T; • Tj) =-\. 
A representative value of 
77 =0.14 (for U = 5.0 and 
Jn = 0.7 eV) is marked by 
dotted line. Parameters of 
the model {21}: / = 35 meV, 
V c = V ab = 0. 




0.20 



S c = (Si ■ Sj) c , S ab = (Si -Sj)ab , (35) 

as well as the orbital correlations, play an important role in the intensity distribution 
in optical spectroscopy. From the form of the above superexchange contributions 
one sees that high-spin excitations (ij) support the FM coupling while the low- 

spin ones, H2 (ij) and H$ '(ij), contribute with AF couplings. 

We have determined the exchange constants in LaVOs by averaging over the 
orbital operators, see Eqs. < fT9b . The case of the ab planes is straightforward as only 
the average densities (n la ) and are needed to determine J i: b, and at large 77 
they follow from the G-AO order in these planes. At 77 = the orbital correlations 
along the c axis result from orbital fluctuations in the ID orbital chain. In this limit 
the orbital correlations are the same as for the AF Heisenberg chain, i.e., (t, • xi) = 
—0.4431 and the ground state is disordered, with (if) = 0. Nevertheless, for this 
disordered state the result for J a t, is similar as for the G-AO phase l25l . 

For the disordered (fluctuating) {a,b} orbital state at 77 = 0, the AF exchange 
interactions in ab planes (see Fig. [3]) result solely from singly occupied c orbitals 
(TT4l) . which are active in these planes and contribute by their double occupancies 
in excited state with AF superexchange. One expects that the exchange constants 
along the c axis in the C-AF phase could be deduced from Eqs. dT9b , as spin and 
orbital order are complementary 1221 . It is quite remarkable that at the same time 
finite FM interactions —J c ~ 3 meV are obtained at 77 = (Fig. |3). They follow 
from the orbital fluctuations which dominate at low values of 77. This mechanism of 
FM exchange adds to the one known in systems with real orbital order at finite 77 
— the latter mechanism gradually takes over when 77 increases and the G-AO order 
develops and reduces the orbital fluctuations. At finite 77 > we used the linear 
orbital-wave theory [12| to determine the intersite orbital correlations (t,- ■ Tj) and 
the order parameter (t?), for more details see Ref. |25l . At 77 = 0.14 representative 
for LaV03, the FM interactions are stronger than from AF ones, \J C \ > J a t,. Indeed, 
this early prediction of the theory Q agrees qualitatively with larger average FM 
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Fig. 4 Kinetic energy Ki, 
(solid lines) for the optical 
subband n and total K\ c > 
(dashed line) obtained from 
the spin-orbital model )2U . 
Filled circles show the effec- 
tive carrier number N^X (in 
the energy range CO < 3 eV) 
for LaVOj , presented in Fig. 
5 of Ref. [SSJ. Dotted line 
shows k[ obtained from the 
MF decoupling l |37> . Param- 
eters: r] = 0.12, V c = 0.97, 
Vat = 0.2/. 
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exchange J c < in the C-AF phase of YVO3 than the AF exchange J a b > in the 
ab planes, see below. 

We emphasize that the strong FM exchange along the c axis follows from the 
orbital fluctuations, and the rigid G-AO order obtained in the limit of strong orbital 
interactions {V a b, V c } ( f20l > would give a much weaker FM interaction, 

Jt AO = -\mJ, 06) 

see Fig. [3] The FM interaction J{?~ AO d36t vanishes at 77 = 0, is triggered by finite 
Hund's exchange T] and increases in lowest order linearly with 77. This behavior 
follows the conventional mechanism of FM interactions induced by finite Hund's 
exchange in the states with AO order, as for instance in KCUF3 J8) or in LaMnOs J2 ■ 
A crucial test of the present theory which demonstrates that orbital fluctuations 
are indeed present in LaV03, concerns the temperature dependence of the low- 
energy (high-spin) spectral weight in optical absorption along the c axis -K[ c) /2J. 
According to experiment (35 j it decreases by about 50% between low temperature 
and T = 300 K. In contrast, the result obtained by averaging the high-spin superex- 
change term H^(ij) ( f29b for polarization along the c axis assuming robust G-AO 
order is, 

^- AO =ln(s c + 2), (37) 

where the spin correlation function s c d35l ) is responsible for the entire tempera- 
ture dependence of the low-energy spectral weight. Equation ( f3Tb predicts decrease 
of w c \ of only about 27%, see Fig. [4] and the maximal possible reduction of 
reached at s c = in the limit of T —> °° is by 33%. This result proves that the scenario 
with frozen G-AO order in LaV03 is excluded by experiment (5). 

In contrast, when a cluster method which allows to include orbital fluctuations 
along the c axis is used to determine the optical spectral weight from the high-spin 
superexchange term ( f29T > 0321 . the temperature dependence resulting from the theory 
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Fig. 5 Spin-wave dispersions 
fflk obtained in the LSW the- 
ory )39t for the C-AF phase 
of YVO3 (lines), and mea- 
sured by neutron scattering 
at T = 85 K GD (circles). 
Parameters: J u t, = 2.6 meV, 
J c = -3.1 meV, S s = 0.35, 
and K- = 0.4 eV (full lines), 
K z — (dashed lines). The 
high symmetry points are: 

r = (0,0,0), m= (w,w,o), 

/S = (ff,w,w),Z=(0,0,w). 




follows the experimental data |35 1. This may be considered as a remarkable success 
of the theory based on the spin-orbital superexchange model derived for the 
perovskites. 

However, the experimental situation in the cubic vanadates is more complex and 
full of puzzles. One is connected with the second magnetic transition in YVO3, as 
we already mentioned in Sec. [2] The magnetic transition at 7/v? = 77 K is particu- 
larly surprising as the staggered moments are approximately parallel to the c axis in 
the G-AF phase, and rotate above 7]v2 to the ab planes in the C-AF phase, with some 
small alternating G-AF component along the c axis |36l . While the orientation of 
spins in C-AF and G-AF phase follow in a straightforward manner from the model, 
i.e., are consistent with the expected anisotropy due to spin-orbit coupling 13T1 , the 
observed magnetization reversal with the weak FM component remains puzzling. 
Therefore, in spite of the suggested mechanism based on the entropy increase in the 
C-AF phase l25l . the lower magnetic transition in YVO3 remains mysterious. Sec- 
ondly, the scale of magnetic excitations is considerably reduced for the C-AF phase 
(by a factor close to two) as compared with the exchange constants deduced from 
magnons measured in the G-AF phase |[T9| . In addition, the magnetic order param- 
eter in the C-AF phase of LaVC>3 is strongly reduced to ~ 1 3jJ.s, which cannot be 
explained by the quantum fluctuations in the C-AF phase (being only 6% for 5=1 
spins [39]). Finally, the C-AF phase of YVO3 is dimerized. Until now, only this last 
feature found a satisfactory explanation in the theory H37II38I , see below. 

We remark that the observed dimerization in the magnon dispersions may be 
seen as a signature of entanglement in excited states which becomes active at finite 
temperature. The microscopic reason of the anisotropy in the exchange constants 
J? c \ = J? c {\ + S s ) and ^f C 2 = ^ c iy — S s ) is the tendency of the orbital chain to 
dimerize, activated by thermal fluctuations in the FM spin chain 1381 which support 
dimerized structure in the orbital sector. As a result one finds alternating stronger « 
^ t -(l + S s ) and weaker °c ^ c {\ — S s ) FM bonds along the c axis in the dimerized C- 
AF phase (with 8 S > 0). The observed spin waves may be explained by the following 
effective spin Hamiltonian for this phase (assuming again that the spin and orbital 
operators may be disentangled which is strictly valid only at T = 0): 
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Fig. 6 Free energies &c (C- 
AF, solid line) and & G (G-AF, 
dashed line) as obtained for 
the spin-orbital model i_ I i us- 
ing the experimental values of 
magnetic exchange constants 
in both phases 1191 . The 
experimental magnetic tran- 
sition temperatures, Tm2 — 77 
K and 7}yi — 116 K, are 
indicated by aiTows. Parame- 
ters: J = 40 meV, 77 = 0.13, 
V a = 0.307, V c = 0.84/. 
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/. L [l + (-l)^]S r S i+1 + ^ £ SrSj+K^istf . (38) 
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Following the linear spin-wave theory 1251 , the magnon dispersion is given by 



OXt(k) = 2- 



2^ + |^ c | + k±^ /2 



1/2 



(39) 



with 



Yk = -(cos k x + cos k y ) , 



cos 2 fe z + 5 S 2 sin 2 £ z 



(40) 
(41) 



For the numerical evaluation of figure we have used the experimental exchange 
interactions |fT9"l : = 2.6 meV, ^ c = —3.1 meV, 8 S = 0.35. Indeed, large gap is 
found between two modes halfway in between the M and R points, and between the 
Z and r points (not shown). Two modes measured by neutron scattering (T9 1 (see 
also figure 1 ) and obtained from the present theory in the unfolded Brillouin zone are 
well reproduced by the dimerized FM exchange couplings in spin Hamiltonian d39b . 
We note that a somewhat different Hamiltonian with more involved interactions 
was introduced in ref. |[T9l , but the essential features seen in the experiment are 
reproduced already by the present model H s with a single ion anisotropy term °= K z . 

As the transition between the two magnetic phases, G-AF and C-AF phase, 
occurs in YVO3 at finite temperature, the entropy has to play an important role. 
As mentioned above, the exchange constants found in the C-AF phase of YVO3 
(Fig. |5]l are considerably lower than the corresponding values in the G-AF phase, 
Jab = Jc — 5.7 meV 1191 . As a result of weaker exchange interactions, the spin en- 
tropy of the C phase will grow faster than that of the G phase, and induce the G^C 
transition. However, starting from our model (fZTJ we do not find this strong reduc- 
tion of energy scale in the C-AF phase. Other mechanism like the fluctuation of % 
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occupancy has been invoked to account for this reduction (25) . Here we will simply 
adopt the experimental values for the exchange constants in the C-AF phase. 

Using linear spin-wave and orbital-wave theory, the spin and orbital entropy 
normalized per one vanadium ion was calculated and compared for both magnetic 
phases of YVO3 l25l . Using the experimental parameters |[T9l one finds that: (;) the 
entropy S?c for the C-AF phase is larger that ^/q for the G-AF phase, and (ii) the 
spin entropy grows significantly faster with temperature than the orbital entropy for 
each phase. Therefore, we conclude that the spin entropy gives here a leading con- 
tribution and is responsible for a fast decrease of the free energy in the C-AF phase 
which is responsible for the observed magnetic transition at Tp/2 E2), see Fig. [6] 



4 Orbital and magnetic transition in the RYO3 perovskites 
4.1 Spin-orbital-lattice coupling 

Experimental studies have shown that the C-AF order is common to the entire fam- 
ily of the RVO3 vanadates, where /?=Lu,- • - ,La. In general the structural (orbital) 
transition occurs first, i.e., 7jvi < Too, except for LaVC>3 with T^i — Too 112011221 . 
When the ionic radius rg decreases, the Neel temperature Tm also decreases, while 
the orbital transition temperature Too increases first, passes through a maximum 
close to YVO3, and decreases afterwards when L.UVO3 is approached. Knowing 
that quantum fluctuations and spin-orbital entanglement play so important role in 
the perovskite vanadates, it is of interest to ask whether the spin-orbital model d2TT > 
is able to describe this variation of Too an d Tjvi with decreasing radius tr of R ions 
in RVO3 ll20l . It is clear that the nonmonotonic dependence of Too on rg cannot be 
reproduced just by the superexchange, as a maximum in Too requires two mech- 
anisms which oppose each other. In fact, the decreasing V-O-V angle (0 along 
the c axis) with decreasing ionic radius rg along the RVOj, perovskites [f40T 443ll re- 
duces somewhat both the hopping t and superexchange J (O, but we shall ignore 
this effect here and concentrate ourselves on the leading dependence on orbital cor- 
relations which are controlled by lattice distortions. 



The model introduced in Ref. [30] to describe the phase diagram of /?V03 in- 
cludes the spin-orbital-lattice coupling by the terms: (/) the superexchange Hj © 
between V 3+ ions in the d 2 configuration with 5=1 spins 0, (ii) intersite orbital 
interactions Hy (120b (which originate from the coupling to the lattice and play an 
important role in the transition between the C-AF and G-AF phase), (Hi) the crystal- 
field splitting oc E z between yz and zx orbitals, and (iv) orbital-lattice term °c gu 
which induces orbital polarization when the lattice strain (distortion) u inreases. 
The Hamiltonian consists thus of several terms ll30l . 



e iR ' Q Tf-gu^ + \NK{u - Mo (#)} 2 • (42) 
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Except for the superexchange fflj ([8]), all the other terms in Eq. ( l42l depend on 
the tilting angle which we use to parameterize the BNO^, perovskites below. It 
is related to the V-O-V angle = n — 2$, which decreases with increasing ionic 
radius r% (0 = 180° corresponds to an ideal perovskite structure). By analyzing 
the structural data of the RVO3 perovskites one arrives at the following 

empirical relation between tr and #: 

r fi = r -asin 2 2#, (43) 

with r = 1 .5 A and a = 0.95 A. 

The crystal-field splitting of {yz,zx} orbitals (E z > 0) alternates in the ab planes 
and is uniform along the c axis, with a modulation vector Q = (tt, 71,0) in cubic 
notation — it supports the C-AO order, and not the observed (weak) C-AO order. 
The orbital interactions induced by the distortions of the VC>6 octahedra and by 
GdFe03 distortions of the lattice, V a b > and V c > 0, also favor the C-AO order 
(like E z > 0). The orbital interaction V c counteracts the orbital superexchange J 
( fTST ). and has only rather weak dependence on so it suffices to choose a constant 
V c = 0.267 to reproduce an almost simultaneous onset of spin and orbital order in 
LaV0 3 , with 7bo ^ T m , as observed ED- One finds r^ p = 147 K taking / = 200 K 
in the present model d42l) . which reproduces well the experimental value 7^j P = 143 
KforLaV0 3 l20l . 

The last two terms in Eq. d42i i descibe the orbital-lattice coupling via the or- 
thorhombic strain u = (b — a) /a, where a and b are the lattice parameters of the 
Pbnm structure, K is the force constant, and is the number of V 3+ ions. Unlike 
E z , the coupling gu > acts as a transverse field in the pseudospin space and fa- 
vors that one of the two linear combinations -^(\ a )i ± \b)i) °f active t2 g orbitals is 
occupied at site i. By minimizing the energy over u, one finds 

geS (0;T) = gu(0;T) = gu (V) + ^(t x )t , (44) 

which shows that the global distortion u(i};T) consists of (i) a pure lattice contri- 
bution mo(#), and (;';') a contribution due the orbital polarization « (t a ) which is 
determined self-consistently. 



4.2 Dependence on lattice distortion 

In order to investigate the phase diagram of the RVO3 perovskites one needs still 
information on the functional dependence of the parameters {E-, V a t,,g e ff} of the 
microscopic model d42l on the tilting angle #. The GdFe03-like distortion is 
parametrized by two angles {#,<p} describing rotations around the b and c cubic 
axes, as explained in Ref. |44|. Here we adopted a representative value of (f> = # /2, 
similar as in the perovskite titanates. Therefore, we used only a single rotation angle 
# in Eq. (l42l . which is related to the ionic size by Eq. (l43l . Functional dependence 
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of the crystal-field splitting E z °= sin 3 $ cos # on the angle # may be derived from the 
point charge model (30), using the structural data for the RVO3 perovskites 140H431 . 
It is expected that the functional dependence of V a t, follows the crystal-field term, 
so we write: 

£-(#) = Jv z sin 3 #cos# , (45) 
Vab (&)=J Vab sin 3 # cos . (46) 

Qualitatively, increasing E z and V a b with increasing lattice distortion and tilting an- 
gle # do favor the orbital order, so the temperature Too is expected to increase. 

A maximum observed in the dependence of Too on rR (or #) may be reproduced 
within the present model d42l i only when a competing orbital polarization interaction 
geg(&',T) (l44t increases faster with # when the ionic radius r# is reduced than 
{E z ,V a t,}. Both mo and (t*) in Eq. (l44l are expected to increase with increasing 
tilting angle Below we present the results obtained with a semiempirical relation, 

geS ($) =J Vg sin 5 # cost?, (47) 

as postulated in Ref. lf30") . Altogether, model (fSTJ depends on three parameters: 
{vzi v ab,v s } which could be selected l30l to reproduce the observed dependence 
of orbital and magnetic transition temperature on the ionic radius r« in the 
perovskites, see below. 



4.3 Evolution of spin and orbital order in RVO3 

Hamiltonian d42l i poses a many-body problem which includes an interplay between 
spin, orbital, and lattice degrees of freedom. A standard approach to investigate the 
onset of spin and orbital order is to use the mean-field (MF) theory with on-site 
order parameters (S z ) (corresponding to C-AF phase) and 



i 1 



(48) 



as well as the coupling between them which modifies the MF equations, similar 
to the situation encountered in the Ashkin-Teller model 1431 . This approach was 
successfully implemented to determine the orbital and magnetic transition tempera- 
ture, Too an d Tjvi in LaMnC>3 J2- It was also applied to the RVOj, perovskites |46] 
to demonstrate that either spin or orbital order may occur first at decreasing tem- 
perature, depending on the amplitude of hopping parameters. However, Such a MF 
approach uses only on-site order parameters and cannot suffice when orbital fluc- 
tuations also contribute, e.g. stabilizing the C-AF phase in LaVC>3 — then it 
becomes essential to determine self-consistently the above on-site order parameters 
together with orbital singlet correlations (l23l on the bonds (ij) \\ c. The simplest ap- 
proach which allows us to determine these correlations is a cluster MF theory for a 
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parameter (t: : )c J48t (dashed 
lines), and spin order pa- 
rameter (S z ) (solid lines) for 
LaV0 3 and SmV0 3 (thin 
and heavy lines). Parameters: 
V c = 0.267, v- = 17, v ub = 22, 
v, = 740. 



Fig. 7 The orbital polar- 
ization (if) (dashed-dotted 
lines), G-type orbital order 
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bond coupled to effective spin and orbital symmetry breaking fields which originate 
from its neighbors in an ordered phase. The respective transition temperatures are 
obtained when (S z ) > ((S z ) = 0) for T < T m (T > T m ), and (x z ) G > ((x z ) G = 0) 
forT <T 00 (T>T 00 ). 

Making a proper selection of the model parameters {v z ,v a b,v g } one is able to 
reproduce the experimental phase diagram for the onset of G-AO and C-AF order 
in the RVO3 family in the entire range of tr, see below. We start with presenting 
an example of the orbital and spin phase transition in LaV03 and in SmV03, see 
Fig. [7] By selecting V c = 0.26/ both G-AO and C-AF order occur simultaneuosly in 
LaVC>3 below Too = Tjvi — 0.737. The crystal field splitting E-, orbital interaction 
V a i,, and the coupling to the lattice g e g are rather small and do not influence the 
order in LaV03. We emphasize that orbital correlations along the c axis are here 
practically as in the AF Heisenberg chain, (t, ■ Xj) ~ —0.44, and the orbital order is 
considerably reduced, (x z ) G — 0.32. The orbital polarization in LaV03 (t x ) ~ 0.03 
is rather weak at Tm , and is further reduced with decreasing T < 7bo- Note that also 
spin order parameter is expected to be reduced below (S z ) = 1, but weak quantum 
fluctuations in the C-AF phase l39l were neglected here. In contrast, in SmV03 the 
phase transitions separate: the orbital transition occurs first at 7qo — 0.86/, and the 
magnetic one follows at a lower T^\ — 0.657. Already in this case the transverse 
orbital polarization is considerably increased, with (t x ) ~ 0.20 at T^\ (see Fig.|7|, 
and further increases with decreasing r% (not shown). Note that the polarization 
(t x ) does not change close to 7bo, and only below T^\ gets weakly reduced due 
to the developing magnetic order, in agreement with experiment J43). The G-OO 
parameter is here stronger as the singlet orbital fluctuations are not so pronounced 
when T -> 0, being (t z ) g ~ 0.37. 

The key features of the present spin-orbital system which drive the observed 
dependence of 7bo and T^i on r# l20l is the evolution of intersite orbital correlations 
are: (i) the gradual increase of the orbital interactions K c ,i,zfz z [Fig.[8ja)], and (if) 
the reduction of orbital fluctuations on the bonds along the c axis, described by the 
bond singlet correlations (t, • Xj) [Fig.[HJb)]. The parameter K a b in Fig.|HJa) consists 
of the superexchange contribution « J ( [T8l and orbital interaction V c ,b < f20t induced 
by the lattice distortion. While the superexchange does not change with decreasing 



Orbital fluctuations in the R VO3 perovskites 2 1 




i-r (A) r R (A) 



Fig. 8 (a) The width of magnon band Wc-af for finite g e ff (circles) and without orbital-strain 
coupling (g e ff = 0, dashed), and orbital interactions in ab planes K u \, (squares) in the C-AF phase 
of cubic vanadates (the points correspond to the SVO3 compounds of Fig. (b) Evolution of 
the orbital order parameter (tf)c (filled circles), transverse orbital polarization (tf) (squares), 
and orbital intersite correlations |{t,- • T ; }| (diamonds) along c axis at T = 0. Parameters: v z = 17, 
Vab = 22, v g = 740. 



rR, the latter term increases and induces the increase of Tqq from LaV03 to YVO3. 
This increase is similar to that observed in the /?Mn03 manganites [;2T|. Thereby 
the bond angle decreases from 157.4° in LaV0 3 to 144.8° in YVO3. 

While the singlet correlations are drastically suppressed from LaV03 towards 
LUVO3, the orbital order parameter (t z )g somewhat increases from LaV03 to 
SmV03 (see also Fig. |7). At the same time the orbital polarization (t x ) increases 
and soon becomes as important as the orbital order parameter, i.e., (t x ) ~ {t z )g- 
Further increase of the orbital polarization towards LUVO3 suppresses the G-AO 
parameter, so (t z )g passes through a maximum and decreases for rR < 1.22 A. 

It is remarkable that the above changes in orbital correlations induced by the 
lattice suppress gradually the magnetic interactions in the C-AF phase, although the 
value of J remains unchanged. This is well visible in the total width of the magnon 
band, Wc-af = 4(J a b + \J C \) (at T = 0) l25l . shown in Fig. ©a), being reduced 
from ~ 1.847 in LaV03 to ~ 1.05/ in YVO3. This large reduction qualitatively 
agrees with the rather small values of the exchange constants in the C-AF phase 
of YVO3 |fl9l , see also Fig. [5] This reduction is caused by the suppression of the 
singlet orbital correlations (t, • T;) by the increasing coupling to the lattice g e ff(#) 
when rg decreases. Note also that this effect would be rather small for g e ff = — 
this behavior is excluded by experiment. 

Following Ref. 03 01 . we argue that the gradual reduction of the orbital singlet 
correlations in favor of increasing orbital polarization is responsible for the evolu- 
tion of the orbital transition temperature Too in the experimental phase diagram of 
Fig. [9] which is reproduced by the theory in the entire range of available rR. The 
transition temperature 7bo changes in a nonmonotonic way, similar to the orbital 
order parameter (t z )g at T = [Fig. Ob)]. After analyzing the changes in the or- 
bital correlations, we see that the physical reasons of the decrease of 7bo for small 
(large) rR are quite different. While the orbital fluctuations dominate and largely 
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Fig. 9 The orbital Too and 
magnetic 7Jvi transition tem- 
perature for varying tr in 
the SVO3 perovskites, ob- 
tained from model d421 for: 
v g = 740 (solid lines) and 
Vg = (dashed lines). Circles 
show the experimentat data of 
Ref. 1201 . The inset shows the 
GdFeC>3-type distortion, with 
the rotation angles # and <p. 
Other parameters as in Fig. 
[8] This figure is reproduced 
from Ref. (30). 
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suppress the orbital order in LaVC>3, the orbital polarization takes over near YVO3 
and competes with G-AO order. 

While the above fast dependence on the tilting angle # of VOg octahedra in 
the family was introduced in order to reproduce the experimentally observed 

dependence of 7bo on rR, see Fig. [9] it may be justified a posteriori. It turns out that 
the dependence of g e ff on the ionic radius r% in Eq. d4Tb follows the actual lattice 
distortion u in RVOt, measured at T = (uq), or just above (ui) DPI . Also the 
orbital polarization (v x ) is approximately oc sin 5 #cos#, and follows the same fast 
dependence of g e ff(# for the RVO3 perovskites (Fig.fTol). This result is somewhat 
unexpected, as information about the actual lattice distortions has not been used 
in constructing the miscroscopic model d2lT) , These results indicate that the bare 
coupling parameters {g,K} are nearly constant and independent of rR, which may 
be treated as a prediction of the theory to be verified by future experiments. 



Fig. 10 Experimental distor- 
tion (in percent) at T = («o, 
circles) and above Tmi (lit, 
triangles) for the RVO3 com- 
pounds 1 40 41 ,43], compared 
with the orbital polarization 
(T t )r=o and with g e ff (47) ; 
geff and g 2 /K are in units of 
/. Squares show the upper 
bound for g 2 /K predicted 
by the theory. Parameters: 
v- = 17, v ab = 22, v g = 740. 
This figure is reproduced from 

Ref. da. 



V 

in 




r„(A) 



Orbital fluctuations in the RVO3 perovskites 

5 Summary and outlook 



23 



Summarizing, spin-orbital superexchange model d42l augmented by orbital- 
lattice couplings provides an explanation of the experimental variation of the 
orbital 7bo and magnetic transition temperatures for the whole class of 
the RVOj, perovskites. A more complete theoretical understanding including 
a description of the second magnetic transition from C-AF to G-AF phase, 
which occurs at Tjs/2 for small ionic radii r« [47 1, remains to be addressed by 
future theory, which should include the spin-orbit relativistic coupling (3T). 

We conclude by mentioning a few open issues and future directions of reasearch 
in the field of perovskite vanadates. Rapid progress of the field of orbital physics 
results mainly from experiment, and is triggered by the synthesis of novel materials. 
Although experiment is ahead of theory in most cases, there are some exceptions. 
One of them was a theoretical prediction of the energy and dispersion of orbital ex- 
citations [12 48 49 1 . Only recently orbital excitations (orbitons) could be observed 
by Raman scattering in the Mott insulators LaTi03 and YTi03 [50, 51]. They were 
also identified in the optical absorption spectra of YVO3 and H0VO3 l52l . The ex- 
change of two orbitals along the c axis in the intermediate C-AF phase was shown 
to contribute to the optical conductivity a {to). 

An interesting question which arises in this context is the carrier propagation in 
a Mott insulator with orbital order. This problem is rather complex as in a spin- 
orbital polaron, created by doping of a single hole, both spin and orbital excitations 
contribute to the hole scattering 11241 . which may even become localized by string 
excitations as in the t-J~ model 1531 . Indeed, the coupling to orbitons increases the 
effective mass of a moving hole in e g systems 1541 . The orbital part of the superex- 
change is classical (compass-like) in ti„ systems, but nevertheless the hole is not 
confined as weak quasiparticle dispersion arises from three-site processes 155115611 . 

As in the doped manganites, also in doped /?i_ T (Sr,Ca) v V03 systems the G-AO 
order gradually disappears 11571 . The C-AF spin order survives, however, in a broad 
range of doping, in contrast to Lai_ x Sr x Mn03, where FM order replaces the A- 
AF phase already at x ~ 0.10, and is accompanied by the e g orbital liquid 11581 at 
higher doping. It is quite remarkable that the complementary G-AF/C-AO order is 
fragile and disappears in Yi_ A Ca A VO3 already at x = 0.02 0571 . The doped holes in 
C-AF/G-AO phase are localized in polaron-like states l59l . so the pure electronic 
model such as the one of Ref. 11551 is too crude to capture both the evolution of the 
spin-orbital order in doped vanadates and the gradual decrease of the energy scale 
for spin-orbital fluctuations. Theretical studies at finite hole concentration are still 
nonexistent in 3D models, but one may expect a transition from a phase with AF 
order to a phase with FM spin polarization at large Hund's coupling, as shown both 
for e g |[T3l and tig ll57l systems. 

A few representative problems related to the properties of RVO3 perovskites 
discussed above demonstrate that the orbital physics is a very rich field, with in- 
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trinsically frustrated interactions and rather exotic ordered or disordered phases, 
with their behavior dominated by quantum fluctuations. While valuable information 
about the electronic structure is obtained from density functional theory [61 1, the 
many-body aspects have to be studied simultaneously using models of correlated 
electrons. The RWOj, perovskites remain an interesting field of research, as it turned 
out that electron-lattice coupling is here not strong enough to suppress (quench) the 
orbital fluctuations l30l . Thus the composite quantum fluctuations described by the 
spin-orbital model (l42l) remain active. Nevertheless, there is significant control of 
the electronic properties due to the electron-lattice coupling. Thus, the lattice distor- 
tions may also influence the onset of magnetic order in systems with active orbital 
degrees of freedom. If they are absent and the lattice is frustrated in addition, a 
very interesting situation arises, with strong tendency towards truly exotic quantum 
states |fl5l . Examples of this behavior were considered recently for the triangular 
lattice, both for e g orbitals in LiNi0 2 EI2 and t 2g orbitals in NaTiO? [63]. None 
of these models could really be solved, but generic tendency towards dimer corre- 
lations with spin singlets on the bonds for particular orbital states has been shown. 
Yet, the question whether novel types of orbital order, such as e.g. nematic order in 
spin models ll64l . could be found in certain situations remains open. 
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